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The relational framework of canonical quantum gravity with non-ultralocal constraints is 
explored. After demonstrating the absence of anomalies, a spatially discretized version 
of the relational framework is introduced. This allows the application of Lieb-Robinson 
bounds to on-shell monotonic gauge-flow when there is a continuous external “time” param¬ 
eter. An explicit Lieb-Robinson bound is derived for the differential on-shell evolution of 
the operator norm of the commutator of discretized Dirac observables, demonstrating how 
a local light cone-like causal structure emerges. Ultralocal constraints do not permit such a 
structure to arise via Lieb-Robinson bounds. Gauge and (3-|-l)-diffeomorphism invariance 
of the light-cone is discussed along with the issues of quantum fluctuations, the nature 
of the nonlocalities, the spatial continuum limit, and possible links to non-commutative 
geometry. 
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1. Introduction and Motivation 


Locality and causality have played pivotal roles in gravitational physics for nearly a 
century. While the concept and important consequences of local light cones are well estab¬ 
lished in general relativity, causality in canonical loop quantum gravity and in spin foam 
models has proven more difficult to elucidate. Almost 20 years ago Smolin [T] pointed out 
the importance of the emergence of classical long range correlations for non-pertubative 
theories of quantum gravity in some appropriate semiclassical limit. These correlations 
are then reflected in how one thinks about causality in quantum field theory on a fixed 
classical background spacetime, specifically that local observables at space-like separation 
commute (micro-causality). During the intervening years, there have been several studies 
of causality in the background independent fully quantum regime, particularly utilizing 
spin foam (covariant) models. However, generally the spin foam models do not support 
causal correlations, unless they are either assumed at the outset mm or involve some kind 
of alteration of the vertex amplitude intended to describe a local orientation [4]. On the 
other hand, approaching the issue of semiclassical causality from the canonical point of 
view is even more conceptually challenging since that approach is a “time-less” formalism. 
Consequently, the consensus expectation is that micro-causality will emerge from some as 
yet to be developed semiclassical limit of quantum gravity, however to date such a causal 
limit for quantum gravity remains lacking. Thus it is mysterious that micro-causality occu¬ 
pies such a foundational place in general relativity, quantum field theory, and the standard 
model of particle physics, yet is still so elusive from a background independent quantum 
gravity point of view. In a broader context, the semiclassical regime of quantum gravity 
is of importance not only from the perspective of causal correlations, but also as a general 
testing ground to examine whether a theory of quantum gravity can behave in familiar 
classical ways in some suitable limit. 

Here we take a simple first step towards understanding these questions starting from an 
unexpected direction. We adopt the relational framework approach to canonical quantum 
gravity which has been developed over several decades 016] HI ID [9]. Then we explore the 
case where the constraints are non-ultralocal, and for reasons to be discussed, limit our 
considerations to the on-shell physics. Next, we apply Lieb-Robinson bounds, originally 
introduced in the 1970s to describe solid-state spin systems [TO], to a spatially discretized 
version of the relational framework and demonstrate how a suitably gauge invariant dif¬ 
ferential local light cone for discretized Dirac observables may be constructed. In essence, 
the local light-cone emerges from on-shell non-ultralocality of constraints in a quite general 
sense within spatially discrete relational framework models with smooth monotonic gauge 
flow described by an external “time” parameter. If the constraints are taken to be ul¬ 
tralocal, then this Lieb-Robinson-based causal structure collapses. Quantum fluctuations 
act to disrupt the local light-cone structure, and a set of general criteria are presented 
which are sufficient for the Lieb-Robinson-based local light-cone structure to survive the 
quantum-classical tug-of-war. 
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The outline of the remainder of the paper is as follows: The relational framework is 
briefly recapitulated in section 2. Section 3 studies non-ultralocality, its freedom from 
anomalies, and introduces “patchy” gauge flow. It concludes with a description of the 
spatially discretized model that is used later on. Section 4 provides an introduction to 
two non-relativistic versions of earlier Lieb-Robinson bounds from the literature: The first 
for Heisenberg operator evolution via time independent Hamiltonians, and another more 
mathematically sophisticated approach for time-dependent Hamiltonians, both on general 
lattices or networks. In section 5 a relativistic differential-time expression for the Lieb- 
Robinson local light-cone is derived for the case of an external time parameter acting as 
the synchronizing conductor of the relational framework’s clock variable symphony or¬ 
chestra. Section 6 discusses gauge invariance and other properties of the relational local 
light-cone. Through a series of questions and answers. Section 7 examines issues related 
to the continuum limit of the spatial discretization, the nature of nonlocality necessary for 
the Lieb-Robinson local light-cone, a possible link to non-commutative geometry, and the 
role of quantum fluctuations. The paper concludes with a brief summary and self-criticism 
of the model in section 8. 


2. Review of the Relational Framework 

In this section we briefly review the necessary points of the relational framework for¬ 
malism. For further details please see mmm- The essential idea behind the relational 
framework is to construct Dirac (gauge invariant) observables from gauge variant (par¬ 
tial) observables. This all starts from the classical phase space ^ description of a re- 
parametrization invariant system whose dynamics is described by a (canonical) Hamilton¬ 
ian that consists entirely of a linear combination of constraints. In the relational framework, 
quantization occurs on the reduced phase space. We start by describing the classical for¬ 
malism. Consider then a set of first-class constraints C/ with I € an arbitrary index 
set. For the case of canonical 4-dimensional general relativity, the index I includes both 
a continuous 3-coordinate index y(I), labeling a point a on the 3-dimensional manifold 
S, as well as a discrete index *(/). The latter index ranges from 0 through W — 1, and 
labels the Nc first-class constraints (gauge conditions) at each point a. These constraints 
satisfy the first-class condition {C/, Cj} = ffj Ck, where generally //j may be a structure 
function, possibly depending on phase space functions. We have assumed all second-class 
constraints have previously been handled by taking ^ to be the surface in phase space 
where they all vanish, and that the Poisson bracket used above is the Dirac bracket. Next 
select a set of gauge variant phase space functions , I G ^ called clock functions or clock 
variables that coordinatize the gauge orbit of any point in phase space within a neighbor¬ 
hood of the (classical) constraint surface (shell) ^ = {m G ^\Ci{m) = 0,V/ G J^}. The 
might include matter (non-gravitational) degrees of freedom. Then Aj = {Cj,T'^} is 
locally non-singular as an (/, J) “matrix,” and one can define the transformed equivalent 
first-class constraints 
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( 1 ) 


C'j = Y,[A-yjCj. 

j 

These obey {Cj,T-^} ss dj, where denotes weak equality, that is equality on-shell, and 
the Hamiltonian vector fields = Xj weakly commute (i.e., commute on-shell). Gauge 
transformations for any phase space function / and set /3^ of reals can be written as 

CXD - 

a^if) = exp(X^) -1 = ^2 • / (2) 

n=0 

x^ = J2p'Xj. 

I 

If one is given a set of real-valued phase space constants (clock parameters) r'^, a weakly 
gauge invariant (Dirac) observable associated with partial (gauge variant) observables / 
and is 

= ^[/](r) = . (3) 

The motivating idea is that ^[/](t) represents the value of / when the clock variables 
take the values r^; i.e., it is a gauge slice or fixing. ^[/](r) is a phase space function, and 
one must compute it first treating /3^ as phase space constants, and only subsequently set 
/3^ = T^ - TG 

One also finds ai 3 {T^) ra + (5^ and ^[r^](r) = on-shell. Also 


on the constraint surface i^[/](t) may be formally expanded as 

One can also derive 

^[/](r) + ^[/'](r) = ^[/ + /'](r) (5) 

mir) mir) ^ ( 6 ) 

{^[/](t),^[/'](t)} ^ « ^[{/,/'}D](r), (7) 

where the Dirac bracket is defined as 

{/, /'}d = {/, /'} - {/, Ci}[A-%{T^, /'} + {/', Ci}[A-%{T\ /}. ( 8 ) 


The formalism simplifies considerably if one can choose canonical coordinates so that 
the clock variables are themselves some canonical coordinates. Then one has a complete 
set of canonical pairs partitioned as {q°',Pa) and {T^,Pi), where the Pj are the canonical 
momenta conjugate to the . Hence in a local neighborhood of the constraint surface one 
can write the constraints as the equivalent set 

Ci = Pi + hi{q%Pa,T-^)-0, 
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( 9 ) 



and setting Pj = —hj{q°',pa,T'^) formally solves the constraints. One can also show that 
the canonical Dirac observables Q^(t) = and Pa('r) = ^[Pa\{T) satisfy equal-r 

canonical Poisson bracket relations. As discussed in mm, the Cj comprise a strongly 
Abelian constraint algebra and obey C'j = Ci + O(C^), hence the Hamiltonian vector 
fields Xj of C'j and Xj of Cj weakly commute. The relations ([5]), Q and ([7|) may be sum¬ 
marized by saying that induces a weak algebra homomorphism w.r.t. pointwise addition 
and multiplication as well as a weak Dirac algebra homomorphism on {/,/^}d- When 
neither / nor f depend on any P/, their Dirac bracket reduces to the Poisson bracket, and 
then there is also a weak Poisson algebra homomorphism. 

Next define 

Hi{T)=Hi{Q%T),Pa{T),T) = ^[h,](r) «/iKQ“(r),Pa(r),r). (10) 

If / is any phase space function depending only on qa,Pa, but not on ,Pj, one has 

-Xfi[f](T)^{H,{T),e[f](r)}- ( 11 ) 

That is, Hi[t) generates the r-parametrized gauge flow of / on the constraint surface. So 
if one specializes to a parametrization invariant dynamical system whose canonical Hamil¬ 
tonian vanishes, one may refer to the Hj{t) as the (r dependent) physical Hamiltonians. 
In the following we will make extensive use of m, and it is important to realize that both 
that PDE as well as its integrability condition hold only on-shell [7]. So we will henceforth 
limit ourselves to on-shell physics. 

In order to quantize the system on the reduced phase space where the classical constraints 
are valid, the gauge invariant canonical variables Pa{T)]T^( t), Pj{t) are mapped to 

operators Pa{T)]T^( t), Pi{t) which generate the quantum algebra if with the usual 

equal-r canonical commutation relations. Given if, its representation (carrier) Hilbert 
space ^ may be generated via the GNS construction employing any positive linear func¬ 
tional (state) on if. Here one is reducing phase space before quantizing; that is, all the 
constraints are satisfied at the classical level. We assume that for all r, the physical Hamil¬ 
tonians Pa{T),T) are represented as densely defined self-adjoint operators on 


3. Non-Ultralocality 

We say hj is ultralocal if it only depends on the canonical fields or their spatial gradients 
(of any finite order) at the point y{I) € U. If all the hj are ultralocal, then {hj,hj} oc 
6{y{I),y{J)) (or its spatial derivatives) and therefore one has {Hj{t), ps 0 for 

y{I) / y{J)- In this case, Hj{t) and Hj{t) have no common q^'^pa^T^ for y{I) 7^ y{J)- 
Ultralocality produces significant mathematical convenience and simplification. Moreover, 
all the commonly used classical gravitational constraint algebras, such as ADM and Holst, 
including possible scalar matter fields, possess ultralocal hj, see dl]. In fact, all the known 
interactions in the standard model of particle physics are ultralocal as well. However, a field 
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theory is not required to be ultralocal, just that measurements so far are consistent with 
ultralocality. Here we keep an open mind, and explore the consequences of non-ultralocal 
hj. One should realize from the outset that the constraint Cj algebra remains Abelian for 
non-ultralocal hj, even though neither hj nor Hj will possess weakly Abelian algebras. We 
will see later on that the quantity for y{I) ^ y{J) plays a crucial role in 

Lieb-Robinson bounds. 

It is sensible to first confirm that the quantum dynamics remains anomaly-free even 
for non-Abelian This means that all the classical gauge symmetries are faithfully 

reproduced in the quantum theory. Han m has proposed the condition that the Hj form 
an Abelian algebra to be used as a definition for freedom from anomalies. It is demonstrated 
below that this is overly restrictive, and we provide an alternative condition for the quantum 
dynamics to be anomaly-free. Following m, one seeks to solve the Schrodinger equation 

( 12 ) 


for a unitary propagator between two Schrodinger states 'I'(t),'I''(t') at initial 

“multi-fingered time” t' and final value r. Let denote the space for r, and let c : M ^ 
be a path from t' to r parametrized by a real-valued “external time” t. We will show that 
is independent of the choice of path c as long as dc^/dt > 0, VI, t. This is a 
mathematical representation of the general covariance (gauge invariance) of the quantum 
dynamics. We will need a little more terminology. The real-valued fields (t) specify one 
T at any given t for each y & and each gauge i = 0... Nc — 1- We will refer to a single 
set of configurations (all y,i) {T^{t)}t for all t as a “slicing.” One configuration (all y,i) 
at one given t is called a “slice” from a slicing. A gauge transformation is then a change 
of slicing from {T{t)}t to {f{t)}t', i.e., a change of multi-fingered time. Independent from 
slicing invariance, t can be smoothly re-parametrized to t' = f{t) (a 1-diff). 

Start by rewriting the r-evolution PDF (|12l) as an integral equation: 


U{t, t') 



dr U{a, t'). 


(13) 


Here u is a variable like r or t' taking a value in ^ whose components are: 


r, for J = I, 
Tj otherwise, 
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0-.J = 


where I is the summation index in (I13p . and r is the (real) variable of integration. Iterating 
this leads to a Dyson expansion: 


J T 7 J Tr ^ t't 


( 2 ), 


I,J •'D 


n=l Do 




1 (n) 

dal 


x{Hj,ia(^^)Hi,{a^^^)---Hj„{a(^^)}, 


(14) 


(15) 


where n is the depth of iteration, and Ii,... G One may concretely picture each 
fixed I as an “/-channel” from t' to r. Each a^^^ is a r-vector (a point in ^), having real 
components a^^\K € J^. The /-subscripts on the Hj{a) correspond to the tj integrated 
over when its corresponding argument aj is integrated. 


Suppose the satisfy for all m and K, where here K ^ plays a role 

like the index I in r/. We refer to this as the ordering hypothesis for the path through 
=;T-space from the initial t' to the final r. Now consider the n-th term in the sum on 
the RHS of (fT^ . Suppose among the n channels /i, ...,/„ G occurs pi times, 

..., Jm G ^ occurs Pm times, 1 < M < n, with 0 < pi,... ,pM < n, and pi -|- • • • +pm = n. 
Under the ordering hypothesis one can recast equation (fTSl) as 



- i//JaW)}. (16) 

Here the channels Ij G ,y,j = 1, • • • ,n are drawn from the (pi x Ji -|- p 2 x <^2 + ■ ■ ■) 
J’s. The sum over the /i,... ,/n has been decomposed into sums over the gauge index 
i{Jk) and integrals over the 3-manifold S with respect to the volume element d^yp(y), 
where y coordinatizes S. The existence of a such a volume element is assured once S is 
taken to be an oriented manifold, and it is not generally necessary to have a metric tensor 
on S for that. However, U{t,t') on the LHS of (|16p must be invariant under smooth 
coordinate reparametrizations of S, and this requires the volume element to be similarly 
invariant for scalar Hj. It is straightforward to see that this can be done if it is possible 
to induce a metric tensor gabiu) on S, where the volume element then takes the standard 
form d^y |det Pafe(y)|^/^. gabiv) does not have to be any physical metric. In fact, if the 
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smooth 3-manifold S smoothly embeds into any Euclidean space E = M”', then the Eu¬ 
clidean metric tensor on E induces a suitable metric tensor on S. The (strong) Whitney 
embedding theorem states that if S is m-dimensional Hausdorff and second countable, then 
S smoothly embeds into E for n = 2m. Taking E to have those properties, one then has 
a volume element on it to render coordinate reparametrization invariant on S, 

moreover this is true for any such invariant volume element. Alternatively, taking E to be 
a paracompact differentiable 3-manifold assures it has a Riemannian structure, whose met¬ 
ric tensor can then be used to construct an invariant 3-volume element in the standard way. 

The notation NCC on the RHS of (jl6l) stands for “non-collisional channels,” and arises 
from the following considerations: Each /-channel has a gauge index i and position y on 
E. i and y are independent degrees of freedom for each channel; however, once the number 
of distinct channels has been fixed to be M, two channels with the same i values can¬ 
not occupy the same position y. That is the channels cannot “collide,” as the associated 
merger or splitting of channels would alter the previously fixed number of channels M, 
so ”NCC” can also be thought of as necessary to avoid double or under counting. The 
{Pk} give the number of r-steps or integrations along each distinct channel Jk- Any single 
^(k) g ^ components: one from the sequence , ■ ■ ■ of real integration 

variables in channel Ji; ...; and any one from the sequence • • • , of integra¬ 

tion variables in channel Jm; and (if not integrated) that component of is set equal to 
the corresponding component of the initial t'. Because of the ordering hypothesis, 
differs from at only one index value J, where then Ij. = J. The 

“r-ordering operator” acts to order the r arguments in each /-channel independently 
so that within each channel J they increase from right to left: ctj"^ < • • • < <7® < (Tj\ VJ. 
Notice that the ordering hypothesis has allowed us to remove the path-dependent limits 
on the multiple r-integrations. 

Equation ()16p sums (averages) over all the paths in ^ from t' to r by advancing mono- 
tonically in r/ within each /-channel independently, stepwise over the n integrations, from 
initial t' to final r, as enacted by the operation. The averaging over the paths from 
the initial to final r-slices is carried out for each fixed n,M, first at fixed (NCC) channel 
configuration {{i{Jk)-,y{Jk))}k and fixed stepping configuration {pk}k by the corresponding 
r-integrations. Then the channel configuration and stepping configuration are separately 
averaged at fixed n, M (the result is independent of which of the latter two averages is per¬ 
formed first). Now suppose we have some “external time” parameter t so that r/ = r/(t). 
As long as the path c\t) = Tj{t) between the fixed initial and final r obeys dc^{t)/dt > 0 
for all / € ^ and relevant values of t, one fulfills the ordering hypothesis. This mild mono¬ 
tonicity condition is consistent with the physical picture of multi-fingered time as “flowing 
forwards everywhere.” The overall sign of the monotonicity condition may be reversed, so 
mutli-fingered time then globally flows ’’backwards,” however the choice of that sign does 


not affect the conclusions. 


By adding over all such ways of channel-wise monotonic advancement from r' to r, the 
overall RHS of (fT6]l is insensitive to a change of slicing (gauge) Ti{t) —)• f/(t) for fixed 
initial t' and final r, and depends only on those initial and final configurations. 

This happens regardless of the commuting or non-commuting properties of the It 

occurs because the monotonicity condition takes care of the required r/-ordering within 
each /-channel separately {Tr becomes a t-ordering), and the operator orderings among 
different /-channels (inside T) on the RHS are averaged over (as a sum over monotonic 
paths or slicings between the fixed initial t' and final r configurations). This absence of 
path or slicing dependence of is the relational framework manifestation of general 

covariance (gauge invariance) for the quantum dynamics: we have obtained freedom from 
anomalies for r-paths monotonic in an external time parameter. Under these conditions 
the propagator U more resembles the familiar one from standard (fixed background geom¬ 
etry) quantum field theory. In the absence of monotonicity, the use of the components of r 
as integration variables, such as in equations (fl^ . (fT4)l . and (fT5]l becomes ill-defined: Some 
/-channels could then have ranges of r/ which are traversed multiple times in both senses 
as t advances. This accords with the intuition that these locally negative lapses among 
generally positive ones are indeed somehow physically anomalous. In fact, a similar notion 
in a different guise was used in [3] to generate “causal” spin-foam vertices and amplitudes. 
Monotonicity may be viewed as the relational framework analog of global hyperbolicity 
on Lorentzian manifolds; however unlike the latter, monotonicity assumes no background 
causal structure. 

If one specializes to the ultralocal case, so Hj{t) and Hj{t) have no common Pa,T^ 
for y{I) 7^ yiJ), and to the gauge diagonal case, for which the same thing occurs when 
i{I) 7^ iiJ), then [Hj{t), Hj{t)] ^ 0 for any / 7/ J. Consequently the Tr-ordered product 
in ()16p factors into separate T^j for each /-channel, and U = ] 3 / Ui, as obtained earlier by 

Han [IT] . 

It is also interesting to compare the T,--ordering in (1161) with the well-known T-ordering 
from standard quantum field theory. Weinberg m gives the following Dyson expansion 
for the S-matrix (operator) in Minkowskii spacetime or special relativity: 

00 „ 

5 = i -I- / d'^Xi ■ ■ ■ d'^XnT{j&{xi) ■ ■ ■ J^{Xn)}, (17) 

n=l 

where H = Hq + V, V{t) = exp{iHQt)Vexp{—iHQt), V{t) = f d^xjS’(x,t) in the in¬ 
teraction picture, with Hq the free-Hamiltonian and V the interaction. This is globally 
Lorentz invariant except for the T-ordered operator product. The T-order of two space- 
time points xi and X 2 (the order of their 0-coordinates) is globally Lorentz invariant unless 
(xi—X 2 )^ > 0 {xi—X 2 space-like), so (flTl) introduces no special frame if (but not only if) the 
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J^{x) commute at space-like distances. While this is often referred to as a kind of causality 
condition, in this sense it arises from the invariance of S that occurs because global Lorentz 
transforms x ^ x' = Ax alter the spatial and temporal components of 4-vectors and so can 
then re-order the T seqnence among the {xn}- By contrast, in the relational framework 
representation just discussed, t and all the ti increase smoothly within both slicings r/(t) 
and fi{t) connecting the fixed initial and final r’s, and the gauge transform is no longer 
simply related to the spatial coordinates y. Thus the relational framework has no built-in 
analog of special relativity’s “causality from global Lorentz invariance.” The Tj- in equation 
(I16p will not generally lead to “causality” by itself, regardless of whether the constraints 
are ultralocal or not. 

The quantity [Hi{T)^Hi{Ty\ for y{I) / y(J) is important in the subsequent sections. 
Therefore we spend some time to study it as well as its classical precursor {Hi{t), H. 

We start by computing classically. 


« ff[{hl,hj}D\{T) = ff[{hi,hj}]{T), (18) 


where the hrst relation has used equation ([7]), and the second equality used that the hi 
do not depend on any Pj inside some phase space neighborhood of the constraint snrface. 
We also have Cj = Pj + hj with {Ci,Cj} = 0 (Recall the Ci form a strongly Abelian 
constraint algebra), and one derives 


0 = {Pi + hi, Pj -I- hj] = {hi, hj] + {Pi, hj) -|- {hi,Pj} = {hi, hj] + 


5hj 

6T^ 


6hi 


(19) 


Thus, 


{Hi{t),Hj{t)] 



Shj 

6Ti 


M- 


( 20 ) 


We use the same sign conventions as Han namely 

{Pa,q''} = (+)' 5 a 

{Pi, hj} = {+)Xpj ■ hj = ^ (21) 

^ Spi 6q^ 6q'‘ Spi 


10 










One also has from (IlSp : 

hj}] (r) = hj}] (r) - hj}] ( 


Consequently, 


E n 

{ki} \ J 

Efn 

Efn 

{kj} \ J 

dHjjr) 

dr^ 


kj\ 

{r^ - 

(r-^ - T-^fj 


kjl 

5T^ 


-hj 

V4 n™"" 1 


('t)- 




6hl 

6T^ 


(t)- 


Combining ([20]) and ([23l) one arrives at 

dHj{T) 

dr^ 


6hi 


{'t)- 


Shj 

6T^ 

Shj 

6T^ 

Shj 

5Ti 


(t) 

(t) 

(t) 


( 22 ) 


(23) 


(24) 


Which canonical variables enter <S/j = Hj{t)}? Since Hi{t) = ^[/i/](r) and 

near shell hj is independent of Prs, so Hj has no Prs and neither does £2ij. How about 
the {q°',Pa) variables? Let / be an arbitrary phase space function solely dependent on the 
(q'“,Pa)’s (and not containing any T^,Pk variables). One then has, 

{Hj{T),ff[f]{T)} « ^[{/i/,/}D](r) = ^[{/i7,/}](r) 

«^[{C7,/}](r)«^[X,./](T). (25) 

In the second equality we have used that hj,f are independent of Pk variables, in the 
third that / contains no variables, and in the fourth that Cj = C'j + 0{C‘^) and 
{fi, 12 } = X■ f 2 - One would like to iterate equation (f25]l to yield 

mir)}} I ff[Xj • /](r)} 

« ff[{Xi . Xj) • /](r) « ffiiXj . Xj) • /](r), (26) 

where in the last equality the weak commutativity of the Xj has been used. If this were 
true, then by using the Jacobi identity and that / is an arbitrary phase space function 
of the q°',Pa variables, one would “find” that Hj{t)} could not be a phase space 

function of any q°',Pa- However, the step marked with ? is invalid because there one is 
using the weak equation (I25j) inside a Poisson bracket. Such a maneuver is inadmissible, as 
no weak equation may be used before evaluating Poisson brackets. Hence {Hi{t),H j(r)} 
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may still be a phase space function of q‘^,Pa as well as the clock variables. 


Suppose the constraints are ultralocal, meaning hj only contains q'^,Pa with y{a) = y{I) 
and clocks with y{K) = y{I). Then 6hj/6T^ oc 6{y{I),y{J)), so by (I2n|) 

0 for y{I) / uiJ)- Hence to obtain {Hj (t), Hj (t)} ^ 0 for y{I) ^ y{J) requires non- 
ultralocal constraints. 


What about spatial parity (P) and time reversal (T) symmetries? Let 
where j = i(/), k = i(J), x = ?/(/), and y = y(J). One has = —.^yi by anti-symmetry 
of the Poisson bracket. One might be concerned that when x ^ y, for ^iy to be non¬ 
vanishing could require violation of P-symmetry (in addition to non-ultralocality); i.e., 
that ^iy could acquire a non-zero spatially odd piece. Such a violation is only necessary 
provided .^iy = J^xy, that is the constraints are “gauge symmetric,” meaning they satisfy 

(Gauge Symmetry) (27) 

which does not generally hold. So a violation of P-symmetry is not generally necessary to 
obtain {P/(r), Pj(r)} ^ 0 for y{I) ^ y{J) if there is non-ultralocality. An examination 
of the formal power series for ^[/](t) dH) similarly shows that ones does not require a vi¬ 
olation of T-reversal symmetry either. The absence of requiring P- and/or T-violation to 
obtain ^ 0 for y{I) ^ y{J) is reassuring since the gravitational interaction 

is not expected to violate those symmetries. 


We now present a very simple toy example of non-ultralocality. Suppose the Hamiltoni¬ 
ans have the special form 

P{x) = h\x)+ [ d^y K^^{q^{x),pa{x),q^{y),q’'{y)]x,y)T^{y), (28) 

Jt, 

with the non-ultralocal term chosen linear in the clock variables for simplicity. We take 
the clock variables as 3-diff scalars and the kernel to be a weight-one 3-density. The 
first term on the RHS is ultralocal and does not contribute to {P/(t), Pj(r)}. One finds 


^(r) = {Pf(r),P^(r)} 


6hf 


5JC 

6T,^ 


(t) 


[K^^{q‘^{x),pa{x),q^{z),pb{z);x,z) - K'^\q''{z),pb{z),q'='{x),paix)-, z,x)]{t). 


(29) 


Next we examine the r and external time gauge flow in more detail. This is the 
bridge we will need to cross to get to the Lieb-Robinson bounds. From the classical 
gauge flow equation (fTT|) and again following Han’s sign conventions that {Pa,Q^} ~ 5^, 
[Pair), Q'^(t)] = (—i){Pa(r), (5^(r)} = {h is set to unity), one has the on-shell 
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quantum gauge flow equation: 

d 


— Ri (i) 


t) 


(30) 


where / may depend on q^,pa canonical variables but not the , Pj types. Adopting the 
Ansatz 

^[/] W = exp [iM(r)] {^[f]{0)) exp [ - fM(r)] (31) 

for some self-adjoint operator M independent of /, one infers from the gauge flow equation 
that 

" dM{T) 


dri 


-Hiir), 


t) 


0 . 


(32) 


Since / is an arbitrary phase space function of q°',Pa, one has that on-shell dM{T)/dr^ — 
Hj[t) may depend on the but not on the q°',Pa or Pk- Here we will assume the simplest 
case, that is 

dM{T) 






which is sufficient, but not necessary. Now decompose M as 

M{t) = 

z 


(33) 

(34) 


where Z is called a “patch” and is just the support of Mz on S. The motivation behind 
this is as follows: Each / € contains a continuous spatial coordinate y{I) as well as a 
discrete gauge index i{I). For nonlocal hj and Hi, besides the canonical variables at y{I) 
(called the “central site”) there are other “nearby” canonical variables living at y' ^ y{I) 
which also enter hj and Hj as their “entourage.” We define the patch Z{I) to consist of the 
central site y{I) together with all those nearby y' where its entourage reside. We assume 
the patches to be bounded and not to take up all of S. Denoting the central site of Z by 
yc{Z) we can set 

Hz = Hz(i) = (35) 

Using (p3]l with fixed /, 

Hz{i){r) ^J2dMx{T)/dT^. (36) 

Later on we will be more interested in the external time t gauge flow of ^[/](t) than in 
the r-flow, so we let 


m = 

I 


dT^{t) 

dt 


Hi{t) ^ Y. 

I 


dr^it) 

dt 


E 


dMxir) 


M{t) = Y^x{r{t)). 

X 


dM{t) 

dt 


(37) 

(38) 
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It then follows from the r gauge flow equation (f30]l that 

A 

dt 




dt 




L I 


dt 


t) 


= {i}mt),mir{t)) 


( 39 ) 


As expected, one sees that H(f) generates gauge flow in external time. More explicitly, 
setting S’[f]{t) = ^[f]iT{t)), ([39]) implies 


liin^[/](t + e) = ^[f]{t) + ie 




+ 0(e^) = exp (ie]HI(t)) S’[f]{t) exp ( —ie]HI(t)). 
We would like to put ]H(t) = YliZ ^z{t) and figure out what Hzit) is, hence by (l37)l 




( dr^t) 
dt 


Hi{t) 


dt ^ ^ dt , - - 

z I z I-.Z{I)=Z z \ / Z{I)=Z 

(41) 

where the last equality follows from the fact that there is only one patch Z[I) with central 
site y{I) corresponding to Z: y{I) = yc{Z). Consequently, 


A^c-l 




Z i=0 




And pulling it all together: 

n{t) = Y,Hz{t) 


z 

ATc-l 




i=0 




(42) 

(43) 

(44) 


M(t) is a “patchy” Hamiltonian generating relational framework gauge flow in external 
time, whose patches Z are based on the non-ultralocality of the original hj. This patchy 
representation of t gauge flow unlocks the door to applying the Lieb-Robinson bounds to 
be introduced in section 4 below. 


Before immediately moving on to the Lieb-Robinson bound, we will need to know a little 
more about [Hxit), Hyit)]. Here we discuss its classical counterpart ^ = {Hx(t), Hyit)} 
to gain some intuition about it before proceeding. 

^ is built from 
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S'*?; 


dhf 

%2 


Qj.y2 


(t) 


(45) 


















Denote 


fll{x,y) = 


dhi 


(t) 


(46) 


We have £2{x,y,t) = = —^{y,x,t) where x = ydX) and y = ydX)-, and 

thus 


^{x,y,t) 


Nc-l 

E 


11,22=0 


dTidx) dTidv) 

dt dt 




(47) 


As a reminder, ^ may be a phase space function of the T^,q°‘,pa variables, but not Pk 
variables, or it could just be a phase space constant, and so include r’s. As a simple 
example, here we will try to construct / from just the r’s. We take /j(x,y) to be {x,y) 
symmetric (P-conserving) and {i,j) anti-symmetric, such as /j(a:,y) = {Ti{x) — Tj{x)) + 
ixi{y) ~ d(i^))- Then for the sum in (HT)) to be non-vanishing, we have to (^ 1 ,^ 2 ) anti¬ 
symmetrize the brackets in that expression to obtain 


^{x,y,t) 


N^-l 

E 

21,22=0 


dTidx) dndy) 

dt dt 


dxidx) dndy) 

dt dt 




(48) 


which is overall (x, y) anti-symmetric as required. So as y —x, ^ 0, but away from 

X = y, ^ IS non-vanishing. We expect it to decay as Hx’s and Hy's patches X and Y 
cease to overlap, but that behavior is not well captured by this toy model for 


Lieb-Robinson bounds were originally intended to study spin systems imbedded in a 
solid-state lattice, so they are naturally discretized. This lattice may be extended to in¬ 
clude a general network and is not limited to a periodic tessellation of 3-space by polyhedra. 
Discretization achieves signihcant mathematical simplihcations, so we will follow that ap¬ 
proach in this initial investigation of relational framework with non-ultralocal constraints. 
We discuss the limitations and issues related to discretization and its continuum limit in 
section 7. A continuum approach will be left for future research. 

Here we describe the discretization of the 3-manifold S into a (generalized) lattice A. 
Associated with each lattice site j € A C S is a D-dimensional Hilbert space. Unless oth¬ 
erwise indicated A will have finite size (cardinality); alternatively, A may be taken to be a 
finite sub-lattice of some countable lattice P C S. The Lieb-Robinson bound does not de¬ 
pend on the dimensionality D, and the Hilbert space for for the entire system is taken to be 
the tensor product of the site-based spaces. Capital Latin letters from the end of the alpha¬ 
bet (previously referring to Hamiltonian patches) will now denote sets of lattice sites, and 
|A| designates the cardinality of X. We say an operator O is supported on a set Y of sites 
if O may be expressed as O = lA\y ® P, where lA\y is the identity operator on sites not in 
Y, and P is an operator dehned on Y. In the following sections we will be most interested 
in the complete (Dirac) observables ^[/]('r), where / is a phase space function containing 
neither nor Pj canonical variables, and use the unitary Hamiltonian patchy gauge flow 
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(I5U|1 . (H5I1 . and ([33]). We assume we can take the discretized Hi{t) Hj{Q°'{T),Pa{T),T) 
as a self-adjoint operator on the Hilbert space which is the tensor product of Hilbert spaces 
over lattice sites j{a) included in its arguments. 


We must take a moment to carefully resolve any potential issues that might arise from 
discretizing operators like the ^[/]('r), and to confirm that the discretized equations behave 
as expected, especially from a gauge-flow point of view. To this end, we have to define what 
precisely is meant by the spatial discretization A of an continuum operator constructed 
as a sum of products of the canonical variables. The discretization map A is defined to 
act linearly with respect to any sum of operators. Acting on a product of continuum 
operators, A annihilates (“apodizes” or cuts off) any product which contains one or more 
factors of canonical variables that are not on the lattice. From this definition follows 
A{AB) = A(A)A(H). The procedure we follow is to ’’Diracify” first by constructing 
the continuum Dirac operator ^[/](t) from /, and then to discretize by acting with A. 
One seeks to demonstrate that the continuum gauge-flow equation (|30ll holds when all 
operators are replaced by their discretized images under A; i.e., that A is a gauge-flow 
homomorphism. This is made easier after one notes that the continuum gauge-flow equation 
([soil is an operator equation with both sides (weakly) equal to the operator corresponding 
to the classical expression 





Ui 


)x,(n(.>fj)‘'') •/. 

j 


(49) 


see [m equation 2.13. Then by restricting the free index I to have y{I) € A, i.e. to be on 
the lattice, and applying the A map, it is straightforward to show that 


d 

dr^ 


A(^[/](r)) 


(i) A{HiiT)),A{ff[f]{r)) 


(50) 


One has to interpret the 1 terms as A(i) = Ia = where Ij is the identity 

operator on the Hilbert space at site j. This way all the J’s appearing in the sums in¬ 
side (|50l) are on lattice, and there is no on-shell operator flow to/from the lattice from/to 
non-lattice-sites. It is also simple to show that [A((7/), A((7j)] ps 0, so the discretized 
constraints are weakly Abelian. Henceforth we drop the A whenever it is clear from the 
context that we are discussing a discretization. 


If are lattice sites on S, the Lieb-Robinson bound require a 3-metric d{i,j). If A,B 
are sets of lattice sites, for future use we define 

d{A, B) = dist{A, B) = min d{i,j) (51) 

diam(A) = max d{i,j). (52) 

i,jeA 

When Lieb-Robinson bounds were first applied to solid-state spin systems, introduc¬ 
ing the static metric d{i,j) was innocuous, however in applying Lieb-Robinson bounds to 
gravitational physics there are several issues of serious concern. Already at the purely 
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classical level, d will acquire a dependence on geometric variables included in the 
so if one has a continuously varying external time parameter t, d{i,j) will inherit a con¬ 
tinuous t dependence as well, while d still describes a discretized 3-geometry of S. We will 
show in Section 5 how this t dependence can be accommodated within the Lieb-Robinson 
bound. Still at the classical level, in the continuum d{x, y) for x,y G T, could be taken 
as the proper geodesic distance between x and y. But once S has been discretized, the 
voxelated classical 3-geometric information and the replacement of the PDE describing a 
geodesic by a finite difference equation will introduce a classical discretization “error” into 
d{i,j). Of course, one expects this classical error to become negligible in the limit where 
the (proper) lattice cell size becomes much smaller than any classical length characterizing 
the 3-geometry. On the quantum level, once the classical phase space functions are mapped 
into operators, the well-known more difficult issues of quantum fluctuations, non-vanishing 
expectations of variances, choices of quantum state, and so on, immediately arise. This is 
most apparent for the quantum clock operators T^: What does it mean for a non-trivial 
operator to “take the value € M”? In the fully developed quantum regime, of course, 
there is not even a well-defined 3-geometry at all, so the best one might hope for is that 
one can find some kind of semiclassical regime or limit that supports or approximates a 
3-metric like d{x,y). For now, we will work at a level (classical or semiclassical) where we 
may safely assume we do have a sufficiently accurate t-dependent d{i,j) on the lattice, and 
discretization error, quantum fluctuations, and semiclassical consistency will be discussed 
later in section 7 after we see what the Lieb-Robinson bound can tell us about relational 
framework operator gauge-flow in external time with non-ultralocal constraints. 


4. Introduction to Lieb-Robinson Bounds 

Here we provide a brief and hopefully self-contained introduction to Lieb-Robinson 
bounds. The definitions and theorems will be presented together with some intuition, but 
we refer the more interested reader to references m and m for the detailed derivations. 

From a pedagogical point of view, it is best to start with the simplest case first: A 
non-relativistic spin system on a 3-D lattice [13]. So consider the 1-dimensional transverse 
Ising model for N spins with Hamiltonian 

N-l N 

H = -jY,SfSf^, + BY,Sf. (53) 

^=1 ^=1 

This spin Hamiltonian has the form H = Hz with Hz supported on Z. Lieb-Robinson 
bounds are most suited to cases where ||i7z|| decays rapidly with diam(Z) > 1. Using (|53]1 
and the metric d{i,j) = \i — j\, we see that the Zeeman term has diameter 0 and the Ising 
(exchange) interaction has diameter one. So ||L7^|| = 0 for diam(Z) > 1, and these are 
examples of “finite range” interactions. There are also other forms of decaying interactions 
such as exponential, and so on. One could also place the spins at the vertices of a graph. 
Then H is again a sum of Hz, each Z being two vertices, with Hz non-vanishing only if 
an edge of the graph links them. In that case the metric d{i,j) could be chosen as the 
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shortest path metric, and gives ||-ffz|| / 0 only if diam(Z) = 0,1. (This is quite different 
from the spin networks usually considered in loop quantum gravity!) 

When discussing these kinds of spin systems it is natural to give operators the (non- 
relativistic) time dependence given by Heisenberg evolution: 

^ = exp[iHt]0’{O)exp[—iHt], (54) 

where for simplicity we have taken H to be (explicitly) time independent. Then one has 
the following m- 

Theorem (L-R): Suppose for all sites i G A one has the L-R condition: 

\\Hx\\ |-A| exp[/idiam(x)] < s, (55) 

X3i 

for some positive real constants Let Ax and By be (bosonic) operators supported on 
sets X,Y, respectively. Then if d{X,Y) > 0, one has 

II [Ajv(t),Hy(0)] II < 2\\Ax\\ ||Ry|| y^exp[-/idist(i,y)] [exp(-2s|t|) - 1] (56) 

iex 

<2||Ay:|| ||.By|| |^| exp[-/idist(X,y)] [exp(-2s|t|) - 1]. (57) 

The physical interpretation of this bound is made especially lucid by an argument due to 
Hastings m, which we reproduce here because of its later importance: Given an operator 
A with support X as above, let B^{X) be the ball radius £ about X: B^{X) = {i ■. 
dist(i,y) < £}. Construct the following operator: 

A‘xit) = j dUUAx{t)U^, (58) 

where one integrates over unitaries U supported on A \ Bi{X) using the Haar measure. 
A^(t) has support Bi{X). Since UAx{t)U^ = Ax{t) + U[Ax(t),U^], one has 

\\A^j,{t)-Ax{t)\\< j dU\\[Ax,U]\\. (59) 

Using Lieb-Robinson bound (f57)l to bound the integrand on the RHS, we see A^(f) is 
exponentially operator norm-close to Ax{t) provided £ is sufficiently large compared to 
2s\t\/fi. That is, a time-evolved operator Ax{t) may be approximated to exponential ac¬ 
curacy by an operator A^(t) supported on Bi{X). Therefore Bi{X) has the interpretation 
of an effective f-dependent support for Ax{t), and the (norm) “leakage” of Ax(t) out of 
the “light-cone” Bi{X) is exponentially small. 

Most commonly the Lieb-Robinson bound is cast into the following form: Suppose the 
L-R condition (1551) holds, then there is a constant vlr that depends on s, ^ such that for 
£ = dist(X,y), and £ > VLRt, 

II [Ax{t),BY{0) ] II < |X| ll^lxll IIBvIl, 
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( 60 ) 



and g{i) decays exponentially with i. From the theorem, vlr = 2s/fi. Ax{t) can be 
approximated by supported on the set of sites within distance I = VLR\t\ of X by 

an error whose norm is bounded by i~^VLR \t\g{£) |X| ||Ax||- Bounds on the leakage of 
information (von Neumann entropy) out of the light-cone were studied in [15]. For Hz of 
finite non-zero range, i.e = 0 for diam(Z) > R for some R, the bound may be further 

improved m- If i? = 1, ll^^zll < J then one finds g{£) decays faster than exponentially, 
roughly g{£) ~ exp(—a£^), for positive constant a. However, if Hz has range 0, the discrete 
equivalent of ultralocality, then g, is undefined since diam(X) = 0, and there is more no 
Lieb-Robinson light-cone. 

The intuition underlying the exponential decay is the following: From the proof m 
one finds that the n-th order term of the exponential comes from a chain Hz^, ■ ■ ■, Hz^ 
such that Zi n X / 0, Zi n Z 2 / 0,..., Zn-i H Zn ^ Zn r\ Y / 0, i.e. a chain of n 
patches Zk,k = 1,... ,n each supporting a local patch of H. Successive H^ are mutually 
non-commuting as their ranges overlap, but more distant ones commute as their supports 
are mutually disjoint. So it is crucial for Lieb-Robinson bounds that [Hz^,Hz 2 \ / 0 for 
Zi n Z 2 / 0 and [Hz^,Hz 2 ] = 0 for Zi D Z 2 = 0 for all Hamiltonian patches Zi,Z 2 . This 
why the (classical) relational framework analog {H/(r), Hj(r)} was studied earlier in sec¬ 
tion 3, where the relation between the Lieb-Robinson bound patchy Hz and the relational 
framework Hj{t) is given by (1441) . These chains of successively overlapping Hamiltonian 
patches generate the effective operator support H^(f) (light-cone). 

While there are clear similarities with some features of relational framework non-ultralocality, 
the relational framework Hamiltonians H/(t) have a non-trivial (and non-unitary) r(f) flow, 
see ()23p . In particular, the simple Heisenberg evolution with a time independent Hamil¬ 
tonian (l54|) does not apply to ^[/](T(t)), and one requires an Lieb-Robinson bound for t 
dependent Hamiltonians. A non-relativistic Lieb-Robinson bound including this possibility 
was derived by Nachtergaele, Vershynina, and Zagrebnov (NVZ) in 2011 [TT], which will 
now be sketched. 

NVZ start with vertices x G F, where F is a countable set of vertices. They assume: 
There exists a non-increasing real-valued function F : [0, 00 ) (0, 00 ) such that 


sup y ^ F[d{x, y)) < 00 , and 
ygr 

(61) 

F{d{x,z))F{d{z,y)) 

sup > - ZZ7-T, - ^ < OO- 

F{d{x,y)) 

(62) 


For ^ > 0 define T)j(x) = exp(—/ix)F(x), so \\F^\\ < ||F||,C^ < C. The Hilbert space of 
states for the subsystem at x G F is JY’x. For finite A C F the Hilbert space associated with 
A is JY’r = The algebra of observables supported on A is 

where is the set of bounded linear operators on J^x- If V A 2 , then identify 

with the sub-algebra (8) 1 a2\Ai of -* 2 ^ 2 ) so C The algebra of local 
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observables is defined as 


(63) 


= U ^A. 

Acr 

The C'*-algebra of quasi-local observables is the norm completion of . The support 
of A € si tv is the minimal set X C A such that A = A! (i) 1a\x for some A! € Ax- 
The generator of the operator dynamics is defined for each finite volume A C T and we 
confine our interest to Hamiltonian interactions (NVZ were also able to include suitable 
dissipative terms). This interaction is such that for each finite X C T and for all t, X) 
is an operator in six and X) = <h(t,X). The evolution map for any finite 

A C T and time t, is a bounded linear map s/t_ —)■ si\, 

^ (i)[-I-(t, Z), A] = >Pz(()(A). (64) 

ZcA ZcA 

The 'l'z(t) are bounded linear maps acting on six, for any X C A such that X D Z, which 
are of the form T^(t) (g) idx\z- The 'l'^(t) have norms that generally depend on X, but 
are uniformly bounded as ||'I' 2 (t)|| < 2||<l>(t, Z)||. Let = ^{C^) be the n x n complex 
matrices. We say a map T G ^(siz) is completely bounded iff, Vn > 1 the linear maps 
T (g) idjUn, defined on siz <g) Mn, are bounded with uniformly bounded norm 

||T||cb = sup ||T (g) idjUnll < oo. (65) 

n>l 

Specifically, HT^Hcb is a map defined on si\,\/A. C T such that Z <Z A, which is independent 
of the choice of A in T. 


To obtain an Lieb-Robinson bound, NVZ make the following two hypotheses: Given 
T, d, F as above, 

(1) For all finite A C T, A't, is norm continuous in t, hence uniformly continuous on 
compact intervals. 


(2) For each A, there exists /j, > 0 such that for all t G M, 


One also finds 



sup sup 


ll^gllcb 

Fi^{d{x,y)) 


< oo. 


I|.%(«)ll < E ll*2WII S E E ll'I'z(*)lkb<Pll6;.|A|l|F|l = M.- 

ZcA x,ycAZ^x,y 

By definition of one has Mg < Mt for s < t. 


( 66 ) 

(67) 


Fix some large time T > 0, and for all A G At^ let A{t) for t G [0, T] be a solution of the 
ODE 

—A(t) = A'a A{t) with A(0) = A. 
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( 68 ) 



Because ||^A(i)|| < Mt < oo, this ODE has a unique solution defined by 7 ^^(A) = A{t) for 
0 < s <t <T, where A{t) is the unique solution of (l 68 t) for t € [s, T] with initial condition 
^(s) = A. We say a linear map 7 : j?/ —for C^-algebras ^ is completely positive 
if the maps 7 O id„ : ^ O ® Mn are positive for all re > 1. Here positive means 

positive algebra elements (i.e. of form A*A ) are mapped to positive algebra elements. 
NVZ showed that the map 'yt,s is a unit preserving, completely positive map. For the one 
parameter group of automorphisms induced by the Hamiltonian generators the 

NVZ version of the Lieb-Robinson bound states: 

There are constants v, /j,, c such that for A € £^x,B G jVy, 


[A,B{t)] II < C(H, H)exp 


n{d{X,Y) -vt) , 


(69) 


where C{A,B) = c||H|| ||i3|| min(|V|,|V|). More specifically, given assumptions (1) and 
(2) above, NVZ’s theorem 2 states that 


[Ais),B{t)] II < (2/C^) PH ||H|| ||E|| min {\X\, |V|)x 

X exp ( - fid{X,Y)) exp (||^'||i,^ - s)) - 1 


(70) 


for X,Y C A,X n y = 0. Notice the bound is uniform over the chosen A. To extend to 
uniformity over T, the definition of ||'I'||i,^ in assumption ( 2 ) above should have the sup 
over X, y altered from A to T. 

In order to adapt this result to relational framework gauge flow, we need to know that 
the derivation introduces the quantity 


CB{X,t) 


sup 


ll=^llcb 


(71) 


where for V C A, SSx is the subspace of ^{jYx) of completely bounded linear maps 
vanishing on the identity. Here NVZ use the cb-norm (in contrast to the standard norm) 
to make the denominator independent of A C T. NVZ’s derivation gives 

CB{X,t)<CB{X,s)+ ril^z(r)||C'ij(Z,r)dr, with (72) 

mm = ll'I'z(*)ll < E F,{d{x,y)). (73) 

x,y^Z 
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One also has that Cb{Z,s) = ||-B|| if Z n Y 7 ^ 0 and otherwise vanishes. Iterating these 
equations produces a Dyson expansion: 

00 - 

(74) 

ni 

n=l 

an = (||'&lk;.)"(C^)"-' F^(d(x,y)). (75) 

x^X,y& 

This implies equations ()69p and ()70p above with the (non-relativistic) spatially uniform 
Lieb-Robinson velocity 

/i // 

for all t G [0,r] (temporally uniform bound). This bound on vlr can be utilized to bound 
the norm-leakage of operators outside the Lieb-Robinson light-cone, analogously to what 
was performed earlier for the time-independent Hamiltonian case. 


5. Lieb-Robinson Bounds for the Discretized Relational Framework 

We now apply the NVZ version of a Lieb-Robinson bound [13] to the relational frame¬ 
work discretized as previously described on some lattice or network. It is essential to handle 
appropriately the fact that the 3-metric d{x,y) is both slicing { t ( s )} and slice s (external 
time) dependent, and the Lieb-Robinson bound should preserve the necessary gauge in¬ 
variance. Strictly speaking, d{x, y) should then be denoted as d{x, y\ { t ( s )}, s) for the slice 
at external time s in slicing { t ( s )}, however we will continue to use the abbreviated form 
d(x, y) for convenience. The reader should bear in mind the suppressed slicing and slice 
dependence. 


The key initial step is to replace NVZ’s Z) by 


Afc-l 


^{t,Z)^Hz{t)= Y, 


dr 


MZ) 


it) 


i=0 


dt 


H: 


jc{Z) 


i^it)), 


(77) 


from the spatially discretized version of (1331) . Equation (ITSp may be iterated as CBiX, t) < 
C^\x,t), where the n = 0 iterate vanishes since X CiY = 0. For simplicity we first 


focus on the n = 1 term: 

di) 




\\'i>zir)\\CBiZ,s)dT 


(78) 


< 2||R||(f-s)^ Y _sup ||iLz(s)||. (79) 

iex Z3j,znYj^& 

The factor 2 on the RHS comes from ||'I' 2 (r)|| < 2 ||L 72 (r)||, where 2 enters from bounding 
||[L7z(r),^]|| by 2 ||L7^(r)|| ||A||. The general idea of the rest of the derivation is to insert 
strategically placed uniforming bounds (sups) after introducing an appropriate factor of 
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Ffj^. The ranges of the sups are also important and have to be selected with care. We also 
use sup{AB) < sup(74) sup{B) for A,B > 0. 


iex,fcGy 


\\Hzis)\\ 

F^{d{j,k)) 


FMAk)) , 


(80) 


and using the RHS of (1771) . 

C^bHxa) <2\\B\\{t-s) ^ ( sup )( sup ) 

j£X,keY Z 3 j,k *2e[0,Afc-l] 

(r(5))|| 




8Tif\i) 


ds 


X (sup)( sup )( sup ) 

y r(s) iiG[0,Wc-l] sG[0,r] 

X ( ( sup )F^{d{j,k)) ], 

\ s&[s,t] j 


I i 

l-'-'n 


FM^k)) 


(81) 


where a bounding sup over slicings {t(s)} has been inserted into the middle factor. Hence¬ 
forth we will take dT^{t)/dt > 0, for all which is just the monotonicity condition 

necessary for freedom from anomalies discussed earlier, so the absolute values in the first 
line may be omitted. We have also “extended” the sup over s in the factor containing ||77|| 
from s € [s,t] to ^ G [OjR]) with 0 < s < t < T, where recall T is some “large” external 
time. 

Next, first bound the first factor, containing the s derivatives of r, by taking an overall 
sup over jc G T, thereby rendering that factor independent of Z. Returning then to the 
middle factor with 1177| |, bring in the sum over Z, and expand j G X, k G V to x,y G A D Z, 
for some chosen A C T of finite cardinality. By expanding J, k to x,y G A, more positive 
terms were added. Hence the sup factor containing ||77 /j(t)|| may be bounded by 


(sup)( sup )( sup )( sup ) V 

t(s) ie[0,Nc-l] SG[0,T] x,yGA z^x,y 


Fy{d{x,y)) 


(82) 


which is now conveniently independent of j, k. Therefore the j, k sum on the far left may 
then be moved all the way to the right to act only on the third factor, containing only 
F^{dij,k)). Thus, 


C'^^^(X,t) < 2 ||R||(t-s)( (sup)( sup )( sup ) 


JcGU ielOjAfc-l] sG[s,t] 


N, 


ds 


Ff(^)(r(s))||' 


X (sup)( sup )( sup )( sup ) 

V r(5) ig[0,Arc-l] se[0,T] x,y&K I^f,[d{x,y)) 


X 


V jex,kGY / 


( 83 ) 
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This may be written more compactly as 


C^^\x,t) <2\\B\\{t-s)Ix,Y{F,) 


( sup ) ( sup ) ( sup ) iV, 
icGF j£[o,7Vc—1] se[s,t] y 


ds 


X \\H\\T,ti , 


where we have set 


Ix,y{F^,) = X] )F^(d{j,k) 

j&Xk&Y 


t{s),s 


Similar to NVZ, we assume there exists a real /r > 0 such that 

\\H\\T,t, = (sup)( sup )( sup )( sup ) 2^ , -^ 

r(s) 2e[o,iVc-i] sg[o,t] Z3x,y Fijt(d{x,y)) 


< oo. 


(84) 

(85) 

( 86 ) 


From ((Ml) we see ||i?||r,/i is a bound temporally uniform over s € [0,T] and spatially uni¬ 
form over the chosen hnite subset A C F. This occurs because A C F is any finite lattice 
containing Z C A, and \\H^k^^\T{s))\\ is independent of the choice of A within F. One 
may extend the dehnition of ||/?||T,/i to be spatially uniform over countable F as in NVZ, 
by changing the sup over x, y in (IMl) from A to F. 


We now examine Ix,Y{Fy) in some more detail. Recall T^(x) = exp(—|Ux)T(x) is a 
positive real-valued, non-increasing function of its non-negative real argument, thus for a 
fixed slicing 


and 


sup F^{d{x,y)) < F^ 

S€[s,t] 


Ix.y(F,) <Y.Y.Fp 

xGX y£Y 


inf fi(x,y)), 
s€[s,t] / 

(87) 

ini d{x,y)). 

sG[s,i] X 

(88) 


Alternatively, in the discretized model under study, X, Y are both hnite sets, so 


Ix,Y{Fy) = sup {d{x,y)) 

xeXyeY 

< sup min (|Ai|, |y|) sup .. E Fy[d{x,y)), (89) 

s€[s,t] y&m(X,Y) ^^m(X,Y) 


where m(X,Y) = Y and M{X,Y) = V if |V| > |y|, otherwise m{X,Y) = X and 
M{X,Y) = Y. Since F)^ is positive, this may be bounded by “expanding” both m{X,Y) 
and M{X,Y) to F to yield 

-^x,y(-F;.) < min (|V|, |y|) exp [ - ^ inf d(V,y)] sup ||F||, (90) 

se[s,t] se[0,T] 

where as a reminder, ||F|| = sup,j,gp T((i(x, y)) < oo contains an implicit s depen¬ 
dence through d{x,y). 
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When one bounds the higher order terms n > 1 in the Dyson expansion (17411 for Cb{X, t), 
at order n one initially inserts n — 1 factors of 


F^{d{x,z))Ff,{d{z,y)) exp [-y{d{x, z) + d{z,y))] 

"S: ^ -r-77-7 t- ^5 


Ff,{d{x,y)) 


exp [ - yd{x,y)] 


(91) 


where the triangle property of the 3-metric on a fixed slice of a fixed slicing has been 
applied. Recall the positive real constant C is defined by (IU2]) . and F^{x) = e~^^F{x). 
There are also n factors of Hz- Expanding CB{X,t) < ||i?|| Yln=i ^n{t—s)'^/n\, one bounds 


fin <(2||R'||r,/.)’^ (sup)( sup ) C'^(fi(r(s),s)) 
^ r(s) sG[0,T] 


n—1 


drtis) 

ds 


( sup ) ( sup ) ( sup ) W 

jcGF iglOjAfc-l] sG[s,t] y 

seb4] j(zx k&Y 

We now assemble all these intermediate steps into the final result. Define 

C;^=(sup)( sup )C^. 
r(s) se[0,T] 


(92) 


(93) 


Because external time re-parameterization (1-diff) invariance will require an Lieb-Robinson 
bound restricted to infinitesimal time increments (f — s) —>■ 0, we set fit = (f — s) —>■ 0 in 
the above expressions. Also define 


( sup )||E(fi(r(s),s)|| < (sup)( sup )(sup) VF(d(x,y;r(s),s)) 
sG[t,t+5t] t{s) se[o,r] x£r 

<sup^E(inf inf (i(x,?/)) = ||E||, 

xer^ H(s)s£[o,t] 
yer 


(94) 


and ||E|| is slicing and slice independent. Then one has the relational framework-Lieb- 
Robinson bound, 


[Axit), By it -I- (5f)] 


< 


Ca 


|R|| ||E||min(|A|,|y|)x 


exp ^ (r(f))^ — exp ^ —/i inf ^d{X,Y'^ 


(95) 


where the on the LHS indicates one is referring to a single slice at some t G [0, T] within 
an arbitrary slicing t(s). As a reminder, X,Y are the supports of respectively, 
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with X n y = 0. We have denoted 





We will refer to the quantity 

VLR = —^ Nc\\H\\t,i_i 

h 

as the relational framework Lieb-Robinson velocity. 


( 96 ) 


(97) 


6. Invariance and Other Properties of the relational framework 

Lieb-Robinson Bound 

The relational framework Lieb-Robinson bound (I95jl has a LHS that refers to the t- 
differential {St) behavior of the norm of an operator commutator between discretized ob¬ 
servables A, B near a single slice of some arbitrary slicing. Typically we take the operators 
A,B to have the relational framework form A(^[/](r)), so they are discretized Dirac ob¬ 
servables. The RHS has many factors, some of which are uniformly bounded over slicings 
and slices, and others that are slicing and/or slice dependent. We now describe those 
dependencies. Recall that a change of slicing is a gauge transformation, and a choice of 
slice s within a slicing is a gauge fixing. By definition (|86p . // ~ l/diam(Z) is chosen so 
that IlLfllT,/^ is finite and includes sups over all slicings r(s) and over all slices at s within 
those slicings. Therefore, in addition to being spatially and temporally uniform, ^ and 
||LI||r,/i are slicing and slice uniform as well. Moreover, = sup.,-^^) sup^gjoj^] is fj. and 
F dependent, and hence slicing and slice uniform too. Thus (2//i) IlLfUr,/^ and then vlr 
are both slicing, slice, and spatio-temporally uniform. The same conclusion holds for the 
prefactor of the exponentials on the RHS of (I95p . 


By comparison with the simple t independent Hamiltonian non-relativistic Lieb-Robinson 
bound (I57p . we see the exponentially damped leakage from a local light-cone is governed 
by the quantity S’ in (fMP above. So when 


5D{t) = inf (d{X,Y;T{s),s 


> Vlr sup 

s£[t,t-\-St] 


OS 


(98) 


there is exponentially small leakage of the operator norm from the local light-cone during 
St. SD{t) is a slicing dependent and slice dependent (slice near t) quantity. It is also t 
reparameterization (I-diff) invariant under t ^ t' = f{t), St —>■ St' = f'{t) St with f'{t) > 0 
for t G [0,r], treating and Hi{t) as (3 -|- 1) scalars, and noting that such a relabelling 
of the slice from t to t' does not affect the 3-metric d on the slice. At a fixed slicing and 
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slice, 6D{t) is also 3-diff invariant because the 3-metric d is, even though the spatial dis¬ 
cretization (lattice) itself is not 3-diff invariant. The factor supg(9r/cls) 6t on the RHS of 
()98l) has the same properties. Hence one has on-shell (3-|-l)-diff invariance of the relational 
framework Lieb-Robinson bound local light-cone. 

As shown by Dittrich mm. on-shell one can classically embed every slicing into a 4- 
manifold with a Lorentzian 4-metric. The tangent bundle of this 4-manifold may also be 
smoothly partitioned in a (3 -|- l)-diff invariant way using that Lorentzian metric to define 
local null directions and so generate a 4-metric based null-cone. The relational framework- 
Lieb-Robinson bound local light-cone should coincide with or bound the 4-metric null-cone, 
but this has not yet been explicitly established. The 4-metric null-cone, however, does not 
address the important issue of observable commutator leakage outside the light cone, which 
is the crux of the relational framework-Lieb-Robinson bound. 

One can also derive a relational framework-Lieb-Robinson bound local light-cone struc¬ 
ture with Vij^snpg{dT /dt)6t replacing vlr supg{dT /dt)6t in a new 6D{t) which is slicing 
independent (containing a suPt-(s)) but is slice dependent (still retaining the sup^gji^^+^^j). 
That is, one fixes some external t parametrization T^{t), and looks at slices within 6t of t 
over all the slicings r(t). But this construction seems less physically natural than the one 
described above, where 5D{t) is both slicing and slice dependent, so we will not discuss it 
further. 

A differential Lieb-Robinson bound local light-cone can be “integrated forwards” in t 
from slice to slice within a single slicing to generate “support tubes” for observables. To 
do this, one constructs operators A^(t, t -\- St) over the external time interval [t, t + St] by 
the Hastings method described in section 4 which are exponentially accurate t-dependent 
supports for a discretized Dirac observable gauge-evolving in external time from t to t + St. 
This then may be iterated for succeeding slices spaced by intervals St ^ 0. This is precisely 
analogous to how one “integrates” null-cones on a curved Lorentzian manifold to generate 
causal curves from the locally Minkowskian geometry. 


7. Physical and Conceptual Questions 

We have explored non-ultralocal constraints with the relational framework and derived 
an external time differential local light-cone structure based on Lieb-Robinson bounds us¬ 
ing a discrete spatial lattice or network model. Several physical and conceptual issues 
about this spring to mind and we discuss those here. This discussion is by necessity less 
mathematically rigorous and considerably speculative in some cases. 

Question 1: Aside from slicing and slice dependent factors like (^dr/ds'^St, are vlr and 
the local light-cone “the same” as the classical spacetime into which the relational frame¬ 
work is embedded varies? That is, basically one is taking sups over a large spacetime to 
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construct VLRSupg{dT/ds)6t. But what happens when the entire spacetime, initial condi¬ 
tions and so on, are altered? Would /r stay the same? If fi changed, then according to (fWl) 
vlr would also change as the spacetime under investigation was altered, a potentially fatal 
physical pathology if we expect vrr to be (or bound) the speed of light. 

One way to avoid such an early demise for the relational framework Lieb-Robinson bound 
would be that l//i ~ diam(Z), the proper “typical size” of single on-shell Hamiltonian 
patches, is a proper length much smaller than the minimum over the classical spacetimes 
of any proper curvature scales Lc{x) they contain. That is, one could interpret 1//U as 
some kind of proper finite range ^ of the on-shell constraints or Hamiltonians. Thus, if 
^ is a microscopic scale relative to classical geometrical scales: ^ -C ra.\B.x Lc{x) uniformly 
for all the classical spacetimes under consideration, then different classical geometries but 
with the same constraints will have the same /i, \ \H\\T,^,Cf^ and thus vrr. This takes the 
constraints to be nonlocal, homogeneous in form, and uniformly bounded; it essentially 
requires a large separation of physical scales, which is common throughout physics. In 
addition, | would also have to be far smaller than any particle physics lengths that have 
be probed so far for there not to have been any evidence yet of nonlocality. But this does 
not mean that ^ has to be on the order of the Planck scale Lp ~ 10“^^ m, but it would 
certainly require ^ << 2 • 10“^® m (1 TeV). ^ would still have to be longer than the scale 
needed to have a well defined 3-metric on each slice, which at least requires that ^ 3> Tp. 

Question 2: What kind of terms in Hj{t) would generate the required non-ultralocality? 
Specifically, would gradient terms (of any order) suffice to produce the non-ultralocality 
for a Lieb-Robinson bound? 

This question is closely tied to discretization. In a typical discretization one replaces 
Vxip{x) —>■ + — where 6 is some kind of lattice constant, and j is a site index 

corresponding to the continuum coordinate x. So Hi containing a point-wise gradient in 
the continuum, such as 'ip{x)\/'ip{x), becomes nonlocal in the discretization since it would 
couple sites j and j + 1- But in the continuum Hj remains firmly ultralocal. Indeed, all the 
well known actions for continuum canonical gravity such as the 3-1-1 ADM decomposition 
or the Holst action are ultralocal, and for those cases ^ = 0, and then there is no more 
Lieb-Robinson local light-cone. 

The answer to the question is no, continuum gradients alone are insufficient for non¬ 
ultralocality. The reason for the negative answer is any discretization perceives a contin¬ 
uum gradient as a lattice constant dependent contribution to ^(6) = l//i(6). An n-th order 
gradient is computed by discretization to give a contribution ^ oc nb, which vanishes in the 
continuum limit 6 —>■ 0. Gradients have no non-zero natural scale in the continuum. Hence 
for a local light-cone to emerge by Lieb-Robinson bounds, the non-ultralocal Hj cannot 
be constructed from products fields and their gradients at a single point. Both ^{b) and 
must be independent of 6 as 6 —)• 0. Physical quantities like vrr and the local 
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light-cone cannot depend on any cutoff scale like b. This is the lesson of the renormaliza¬ 
tion group for background dependent quantum held theory, and also requires discretization 
independence in the continuum limit for gravitational theories. In fact, as found in m, 
discretization independence implies nonlocality in 4D discrete quantum gravity. 

This question and its answer lead us to ask ... 

Question 3: What is the role of the discretization in the Lieb-Robinson bound? 

Contact with physical reality occurs when 6 —>■ 0 because any non-zero b spatial dis¬ 
cretization by itself breaks 3-diff invariance since the lattice or network is not generally 
mapped into itself by a 3-diff m- There are, however, perfect actions for discretizations 
that do recover the requisite invariance as 6 ^ 0 for several models [l8]. Even for these 
special cases, the discretization will coarse-grain or voxelate the metric information at least 
over the scale b. This classical error makes the discretized d{j, k) acquire a b dependence. 
Again we have to assume that hmb_>.o(i(x,y; 6) is well-dehned so VLRS\x\ig{dT/ds)5t and 
5D{t) are all also well-dehned in that continuum limit, where the classical voxelation error 
disappears. 

Question 4: How might the required non-ultralocality of hj, Hi{t) arise? That is, if one 
starts from some bare classical constraints specihed by Hf that are ultralocal, how can 
one end up with effective quantum constraints generated by a non-ultralocal 

The simple answer is that the real quantitative origins of non-ultralocality hi he beyond 
the scope of the present work. Indeed, this is like asking what atomic physics lies behind 
the exchange coupling constant J in the Ising model (f5^ . The quantum algebra it at the 
kinematic level is generated by Qa{T),Pa{T) through their canonical commutation rela¬ 
tions. From that algebra we can hnd its representation carrier Hilbert space by the GNS 
construction from any positive linear functional (state) p on il. The quantum constraints 
are p[Pi + hf^] = 0, i.e. the quantum fluctuations satisfy the effective constraints in the 
mean (as an expectation value). A similar process occurs in background-dependent quan¬ 
tum held theory where quantum corrections due to loops are taken into account using a 
held dependent effective (“dressed”) action r[4>] rather than the bare classical action /[$] 
[19]. Connected vacuum-vacuum quantum held theory amplitudes can be computed using 
tree-level (mean held) Feynman graphs with vertices using r[<h] instead of /[4>]. A sim¬ 
ilar “dressing” by quantum huctuations might render non-ultralocal. Unfortunately 
the quantitative details are out of present computational reach for background independent 
quantum gravity without the well-dehned path integral technology of quantum held theory. 

A more speculative answer is that we may not know the “true” physical action or con¬ 
straints for gravity at very short but not yet Planckian lengths, only that they look local 
as far we can tell from our experience at long length scales (above 2 • 10“^® m or energies 
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up to 1 TeV). In that case, the Lieb-Robinson local light-cone would be a long length scale 
manifestation of non-ultralocality of those otherwise inaccessible short length scale con¬ 
straints, a hint that we might not be aware of some deeper physics. One possibility is that 
non-commutative products such as the Moyal-Weyl-Groenewold ^-product and deformed 
diffeomorphisms could play a role. While non-commutative field theories still have La- 
grangian densities and Noether currents, the products of objects in the Lagrangian density 
are nonlocal [20]. The scale of that nonlocality might be associated with ^ above. However, 
non-commutative geometry is not a quantization of the underlying manifold in the sense 
that it does not promote phase space functions to operators. So non-commutative mani¬ 
folds do not describe quantum fluctuations of geometry, and they are on-shell descriptions. 
Instead they require that the geometric and matter field actions be invariant under de¬ 
formed diffeomorphisms of the non-commutative manifold, which are nearly the standard 
diffeomorphisms normally used to describe manifolds like S that leave the usual local action 
invariant. At length scales larger than the nonlocal effects, the noncommutative symme¬ 
tries and constraints should approach the standard ones while an on-shell Lieb-Robinson 
light cone emerges from a relational framework based on the non-commutative action. It 
is worth noting in this regard that the nonlocality induced by non-commutative geometry 
is not of the gradient (spurious) type. 

Question 5: In the relational framework-Lieb-Robinson bound which variables are quan¬ 
tum and which are classical? 

The operators A, 13 appearing in the Lieb-Robinson bound analysis include those con¬ 
structed from Qair), Pair), which may describe either geometrical or possibly matter de¬ 
grees of freedom. These are treated fully quantum mechanically from the on-shell gauge- 
flow (HH). The Tj are real parameters and are never promoted to quantum operators, so 
they do not acquire fluctuations. But the clock variables are promoted to operators T^, 
making the classical relational framework notion of ’’T^ taking the value r/” more problem¬ 
atic. Nevertheless, these clock operators lie buried inside the Hamiltonians Hj, which are 
then treated fully quantum mechanically. More crucially, at the heart of the Lieb-Robinson 
bound lurks d{x,y), the on-slice 3-metric, computed from geometrical variables. In d{x,y) 
those variables, such as the ADM 3-metric Qab, are treated classically, while in A{t),B{t) 
those same variables are treated quantum mechanically. So there seems to be an inconsis¬ 
tency. That is, inside A[t) and B{t) gab is treated as an element of a non-commutative 
(quantum) C'*-algebra, but when ijab enters a computation of d{x,y,T{s), s) it treated as 
an element of a commutative (classical) C'*-algebra. However if we use classical geodesics 
as a basis for computing d{x,y) as an extremum over the quantum expectation (gab), then 
the Lieb-Robinson local light-cone quantity 6D{t) fsee lMj) is mean field (expectation value 
level) with respect to gab over Within the Hamiltonians gab is treated as a 

quantum operator. Thus we can say the Lieb-Robinson local light-cone construction is at 
least mean field with respect to fluctuations entering gab, and it must be stabilized against 
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those fluctuations if it is to survive. 


Question 6; What parameters, if any, delineate a window of survival for the mean-field 
Lieb-Robinson local light-cone in its precarious perch among the tensions between the clas¬ 
sical and quantum worlds? 

To answer this, we introduce 5d as the largest rms quantum fluctuation of the the proper 
lengths d{i,j), and also take the classical discretization error in proper lengths to be be 
roughly the same as the lattice constant b. Lc denotes the shortest classical curvatnre scale 
of any 3-geometry nnder consideration. We also have ^ = 1/^, the range (proper patch 
diameter) of the non-ultralocal constraints. Since a classical metric d{i,j) loses physical 
meaning when d{i,j) < 6d, and quantum fluctuations could play a role in the origin of 

< Q We can then qualitatively delineate fonr distinct physical regimes: 

(a) (5d < ^ < 6 <C Tc Here the discretization is too coarse to resolve the finite range ^ of 
the Hi, which then appears to be ultralocal. The Lieb-Robinson local light-cone does not 
emerge. 

(h) 5d < b < ^ Lc Here the Lieb-Robinson local light-cone is stabilized against the 
quantum fluctuations 6d, the discretization can resolve the non-ultralocal range ^ of the 
Hi, but it does not resolve the quantum fluctuations. The local light-cone is b (discretiza¬ 
tion) independent once 6 <C Q 

{c) b < 5d < ^ Lc The proper range ^ is still separated and immune from qnantum 
fluctuations, while the local light-cone is still discretization (6) independent. 

(d) 5 < Jd ~ ^ -C Lc Now ^ is quantum-limited, a fluctuation limited patch size. We still 
have discretization independence once b < 5d, where the discretization error is no longer 
physically relevant. 

In cases (b) and (c) the Lieb-Robinson local light-cone will survive provided: (1) the con¬ 
tinuum limit 6 ^ 0 is well-behaved, and (2) 5d <C 6D(t) (see [98]). This means 5d does not 
significantly affect the width of the local light-cone at t,t + St, and thereby limits St from 
below. That is, the t-dependence of the classical 3-metric d{x, y) cannot vary too quickly 
in external time, so one can define a mean field or classical differential local light-cone, (b) 
or (c) could correspond to the non-commutative geometry scenario for nonlocality over the 
scale ^ in the presence of a metric solution. In the marginal case (d), where the scale of 
nonlocality is about that of the quantum fluctuations, the survival of the Lieb-Robinson 
local light-cone is too close to call. Unfortunately, we really do not know which regime 
we live in. Bnt if non-commutative geometry provided an action invariant nnder deformed 
diffeomorphisms and nonlocal on scale then that would be on-shell, and could naturally 
separate ^ from quantum length fluctuations Sd. 
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An alternative approach to relieve the classical vs. quantum tension inherent in d{i,j) 
is to use semiclassical (coherent) states 'h [21] of 3-space. The idea is that each length 
on the initial slice t has a quantum expectation {d)\ii = diii, a quantum fluctuation {{d — 
and the classical length del > 0- It is also possible to use other geometric 
quantities besides d{x, y), such as the areas of triangles or volumes of tetrahedra. To achieve 
semiclassical consistency for d{i,j) in a state 'h, one would require for all sites i,j 


[del — 
del 


<C 1 and 


o'd(T) 

del 


<C 1. 


(99) 

( 100 ) 


This makes the notion of classical distance insensitive to the quantum fluctuations from 
the state T. Such a construction would also encounter difficulties in case (d) above, where 
the range ^ is quantum limited. 


8. Summary, Self-Criticism, and Conclusion 

In this work we have explored the consequences of non-ultralocal constraints within the 
context of the relational framework of canonical gravity. It was shown that this leads to an 
on-shell non-Abelian algebra for the physical Hamiltonians, while the constraint algebra 
remains Abelian. Unitary propagators stay anomaly-free for smooth monotonic gauge flow 
in an external time parameter t. A set of Hamiltonians that generate operator gauge-flow in 
t with finite-ranged support patches was derived. After introducing a spatial discretization, 
Lieb-Robinson bounds were reviewed and applied to demonstrate an on-shell differential 
time local light-cone. This local light-cone has the properties that there is exponentially 
small norm leakage of discretized Dirac operator commutators outside the local light-cone, 
it displays suitable gauge (slicing) and (3-|-l)-diffeomorphism invariance, and the local 
light-cone can be “integrated” into “support tubes” for discretized Dirac operators that 
resemble familiar causal curves from general relativity. This entire Lieb-Robinson bound 
local light-cone structure collapses for ultralocal constraints. Therefore non-ultralocality 
together with Lieb-Robinson bounds go an unexpectedly long way towards explaining how 
the standard quantum field theory version of micro-causality, where local observables com¬ 
mute at space-like distances, emerges from the (semi-) classical relational formulation of 
canonical gravity at length scales greater than that characterizing the nonlocality. 

Within the application of quantum field theory to fixed curved background spacetimes, 
one can derive the familiar causal advanced and retarded propagators as inverse wave op¬ 
erators (Green’s functions) for matter fields such as scalar bosons and so on. These show 
that the vacuum expectation values of commutators of canonical fields vanish outside the 
past or future light cone. However, while this is straightforward for Minkowski spacetime, 
to obtain unique solutions for general curved space times one imposes the stringent require¬ 
ment that the manifold be globally hyperbolic. The Lieb-Robinson approach, on the other 
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hand, requires no such corresponding ab initio strong global background causal structure 
assumption. 


The criticism of the Lieb-Robinson bound route from non-ultralocality to local light- 
cones is abundantly clear from the responses to the questions in the previous section. 
While some issues, such as how two different classical spacetimes can share the same Lieb- 
Robinson velocity, or the role of field gradients, are quite clear, many deeper concerns 
remain only partially clarified, or just display our glaring ignorance. These harder nuts 
to crack include; What are the origins of nonlocality? What is the detailed microscopic 
meaning of the range or correlation length ^ = 1/^? What if ^ is about the size of quantum 
fluctuations (marginal case)? What are the specifics of the semiclassical limit or choice of 
quantum states necessary to ensure that quantum fluctuations do not destroy the local 
light-cone? Is it possible to handle the continuum limit more thoroughly than simply to 
assume the required limit is well-behaved? Might non-commutative geometry or field the¬ 
ory play a role in these issues? Each of these questions challenges us to probe more deeply 
into the “atomic” theory underlying the model of condensed matter ancestry presented 
here and stands as motivation for future work. 

Nevertheless, it remains surprising that aspects of causality may be linked to non- 
ultralocality. Adopting ultralocality uncritically might be somewhat like what occurred 
in the 1950’s with parity; A beautiful symmetry, but Nature could be a lot more interest¬ 
ing if She broke it once in a while. 
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